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ON CONNECTION BETWEEN VALUES OF RIEMANN ZETA 
FUNCTION AT INTEGERS AND GENERALIZED HARMONIG 

NUMBERS 

PAWEL J. SZABLOWSKI 


Abstract. Using Euler transformation of series we relate values of Hurwitz 
zeta function at integer and rational values of arguments to certain rapidly 
converging series where some generalized harmonic numbers appear. The form 
of these generalized harmonic numbers carries information about the values 
of the arguments of Hurwitz function. In particular we prove: V/c € N : 

nk —1 ' fit) 

= 2 ^- 1 —f —’ where Hn are defined below generalized 

harmonic numbers. Further we find generating function of the numbers C(^) = 


1. Introduction 

First let us recall basic notions and definitions that we will work with. By the 
Hurwitz function C(s,a) we will mean: 


C(5 




1 




considered for Res > 1, Rea S (0,1]- Function C(S;1) is called Riemann zeta 
function. We will denote it also by C,{s) if it will not cause misunderstanding. It 
turns out that both these functions can be extended to holomorphic functions of s 
on the whole complex plane except s = 1 where a single pole exists. Of great help 
in doing so is the formula 


(l.I) 


C(s) 


2S-1 
2 '*-! - 1 


OO 


E 


j" 


that enables to extend Riemann zeta function to the whole half plane Res > 0. 
We will consider numbers: 



f- (-ly 

^ {mj + iY ’ 


for m S N and i G {I,..., m — 1} . Notice that 1)-^ ^and 

= 7r/4. The number (27+i)^ called Catalan constant K. 
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It is elementary to notice that 

V(2to)) - C(fc, 1/2 + i/{2m)). 

The main idea of this paper is to apply the so called Euler transformation that 
was nicely recalled in [5]. As pointed out there we have: 


k—1 n=0 


where {cLk\k>i ^ sequence of complex numbers and the sequence A"afe is defined 
recursively: A°afe = A^Ofc = A”-iafe - A”-iafe+i = 

Sondow in presented general idea of applying Euler transformation to Riemann 
function. He however stopped half way in the sense that he calculated finite dif¬ 
ferences A" applied to (j -|- 1)“'* only for s being negative integers. We are going 
to make a few steps further and calculate these differences pointing out the role of 
the generalized harmonic numbers in those calculations. 

The paper is organized as follows. In the next section [5] we present an auxiliary 
result that enables application of Euler transformation to the analyzed series. Fur¬ 
ther we present transformed series approximating numbers In Section[3]we 

calculate generating functions of certain series of numbers and functions. More pre¬ 
cisely we calculate generating functions of the generalized harmonic numbers that 
we have defined in the previous section. We also calculate generating function of 
the series of the generating functions that were defined previously. It turns out that 
this calculation enables to obtain the generating function of the series sums that 
appear on the right hand side of (ini. Finally in the last Section 0] are collected 
cases when exact values of numbers M are known. 


2. Euler transformation 
To proceed further we need the following result. 


Proposition 1. Let us denote = Ej=o(-l)^ (”)/(”^/ + R = 0,1, ■ ■ •, 


and the family of sequences defined recursively: = I, Ho,fc = k > 1, 


Vn,fc>0:Ri7)=E-= 


B 


{m,i) 




We have then: 


where 


V m e N ; = 1, = 0, 

U.U ’ n,U ’ n,i m( 2 /m)n+i ' 

{a)n = a(a + 1)... (a + n — 1) is the so called ’rising factorial’. 

V n > 0, fc > I we get: 


.(m.i) 

^n,k 


n! 

m{i/m)n+i 
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Proof, i) The fact that An,o = 0 follows immediately properties of binomial coeffi¬ 
cients. Notice that we have 


^n+l,k 


min ■ 


1 ) 


m(n -|- 1) -h i 
m(n + 1) 


3=0 ^ ^ ^ 


(- 1 ) 


n +1 


Tfl 


(n + -1-1) -I- i 


i=o 

(_l)»+i 
{m{n -I- 1) -I- i)^ 


3=0 

n+1 

3 


(m{n -I- 1) -I- z) 

(_l)n+l 




(m(n -I- 1) -I- z) 

min+iy. 


1 


n +1 


z(rz -I- 1) 




rz -I- 1 
3 


j\{n — j)\(rn{n -|- 1) -|- z) 
1 

11 -I- z 

3=0 


m(n -I- 1) -I- j 


rz -I- 1 
3 


-)/(mj + i)^ 
/(mj+zf-^ 


1 


ta: 


(m,f) 

n+l,k—l 


V . , + 1) 

-mrr n rn(n+l-+ a _ jm+i fW nntinr Ibnl i-r hnvr m(n+l) .(m.i) 

Since ti - ) - m+mn+r ' notice tnat we nave ^(„+i)++„.i 

= 0 from which immediately follows that = ra{i/m ,) —f since = i. Now 

divide both sides of the identity A+;>^ - = tSTi 

Ai~>, and denote b!+I = Ai+'/A+'.We get = 


_ A{m,i) 1 

mn+m+2 ^n+l,/c —1 


B. 


{m,i) 


• 4 ( 771 .i) 

since ++ 1)1 


m(n+l) 4(771,i) 


= l/i 


mn+m+2 n,^ 


. Hence B 


(77l,i) 


E 71 

.7 = 


B) 


(m,z) 


j—0 7nj-\-i'^j,k—l 


k-1 


m(n+l)+i n+l,fc —1 

since Vfc > 1 : 

□ 


Remark 1. In the literature (compare e.g. [I], [2], [l]j there function notions of 
harmonic and generalized harmonic numbers defined by 1/i^t rz > 1. 

Numbers hn '^ are called simply (ordinary) harmonic numbers. 


We are going to dehne differently generalized harmonic numbers. 

Definition 1. For every k gN numbers < Hn^'^ \ defined recursively by Hn'^ 

1 J n>l,k>0 

= 1, ^Vjj zz > 1 will he called generalized harmonic numbers of 

order k. 

Remark 2. It is easy to see that and that Hn '^ is an ordinary n—th 

harmonic number. 

(k) 

Remark 3. Notice that Hf ^ is a symmetric function of order k of the numbers 
{1,1/2,..., 1/rz} hence it can be expressed as a linear combination of some other 
symmetric functions of order less or egual k. For example we have: Hn^ = hn'^ 
= Hn (the ordinary harmonic number), Hn'^ = Ilf/2 + hn'^ (2, Hn^ = n))/6 -I- 
Hnhn'^(2 + hn'^/3 and so on. 


Remark 4. Notice also that recursive equation that was obtained in the proof of 
Proposition Q] z. e. 

_ m(n + 1) _ _1_ ■ 

m{n + \) + i^ ~ m{n + l) + i 
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is valid also for fc = 0, —1, —2,. Of course then we apply it in the following 

form: 

An+i,k-i = {rn{n + 1) + i)An+i^k - m{n + 

getting for example : A^^di = 1, A^d^i = —to, = 0, Aq "^2 = = 

—m{m + i + 1), A^^_f ^2 = 2to^ , ^^-2 = 0 n = 3,4,... . The fact that 
= 0 for n > k + 1 was already noticed, justified and applied by Sondow in [S] . 

As a corollary we have the following result: 


Theorem 1. 


( 2 . 1 ) 


■4) 


U<J ■ 

^ 2"+lTO(l/TO)n+i 


B, 


{m,i) 


where numbers are defined above, 

i) In particular : 


( 2 . 2 ) 


hk-\-l 


^ 1^(2,1) _ 2fc+l (-l)^^2fc 

to 2^+1 2fe+l 2(2TO)2^+l(2fc)! ’ 


(2.3) 


f'-> 


= k=£ 


n—0 


n!(H2n+l-gn/2) 

2(2n+l)!! 


where i?„ denotes n — th (ordinary) harmonic number, 
a) for m = i = 1, k € N: 


(2.4) 


E 

1=1 


(-1)^-1 


= E 


H, 


(fe-i) 


n2" ’ 


and eonseguently for fc = 2, 3,... 


(2.5) 


■yk—l 


c(fc) = . E 


H, 


(fc-i) 


2fe-i - 1 41^ n2" 

n=l 


Proof. Applying Euler transformation to the series we have 

00 

n—0 

Now it remains to apply Proposition [I] i) To see that (12.11) reduces to (12.31) when 
fc = 2 , TO = 2 and i = 1 notice that i3„ g = 1 and consequently = 

I]"=o l/(2j + l) =i?2n+i-2n„. Further we have (l/2)„+i =n”=oO + l/2) = (2»^+ 
l)!!/2”+^. To justify (12.21) we have to observe that 2M2^’)|)’™^ = S(2k + 1, 2m, to) = 

(-1)^ , _ 1 (-1)^, The fart that (- 1 )'’ 

l^j=-ao (j2m+m)2'=+i “ + l (j2+l)2'"+i ' ^ Z^j=-oo (j2+l)2'=+i 

_ back to Euler and was recalled in [^. 

ii) If TO = i = 1 we have (!)«+! = (n + 1)!. Recall also that then = H^+i- 

(12.51) follows additionally (II.IE □ 
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Remark 5. Notice that when i = \ then the sequence generated by the 

reeursion: = 1, + !)• arguing by induction 

we see that Vn > 0; Consequently we deduce that the sequenee 

is increasing which is not so obvious when eonsidering only definition 
of these numbers. It is also elementary to notice that 

lim = 1. 

k — 

In partieular we deduce that the sequence {C(A:)(1 — l/2^“^)}fc>i is increasing. 


Remark 6. Notice that one ean easily prove (by induction) that \/n,k G N : 1 < 
< n. Hence, utilizing 1^.51) we have: 


ln2-l/2 
2™+i(m + 1) “ 


m 


ofc-l ™ rr(fc-l) 
_f_ ^n+1 

2 fc-i - 1 ^ 2"+i(n + 1) 

n—Q ^ ^ 


< 


1 

2m+l ’ 




since 2 ^- 1 -i — 2 for k <2 and further 

1 / 2 "+^ < and ^ 

E oo 1 _ ln2—1/2 

n—m+1 2»^-'^+i(n—m+1) 2"^+^(m+l) ' 


rrik — 1) 
^n + 1 

0 2^+i(n+l) 


< 


m+l 


> 


ondEZ- 


n=m+l 2'‘+i(n+l) 


— 2'"+i(m+l) 


Remark 7. Formulae i2.1\) and i2.5\) can be considered as a series transformation 
to speed up its convergence. Apery for C(3) in his breakthrough paper and later 
Hessami Pilehrood et al. in [3] obtained series transformations to speedup series 
appearing in the definitions of Riemann or Hurwitz zeta functions. As it is remarked 
in [3] all these transformation give series more or less of the form Cnjwhere Cn 
= 0(1), but for different arguments of f one gets very different series in a very 
different, particular way. Apery’s one is one of the simplest. Formulae mi and 
ifO) offer unified form of the transformed series and speed of convergence is only 
slightly worse. Namely of the form c„/2". 


Remark 8. Notice also that analyzing the proof of Proposition\^we can formulate 
the following observation. Let us denote A^Z = YZj=oZZ iZ/Z) + Z 
Re(s) > 0. Then 

Am,:) _ rn{n + 1) Am.i) _ __ Am,:) 

m{n + l)Fl ~ m{n+l) + l 

Hence keeping in mind that A^Z ~ 1/^* assuming that we know numbers 
I ZZ-i f we are able to get numbers | A^Z )• and consequently find C,{s,l/m) 
In particular if m = I = 1 we get An+i,s — = ;;/^^n+i,s-i where we 

denoted An,s = aZP to simplify notation. Consequently we deduce that An^a = 
ZiZiZi Sinee we can iterate this relationship we see that the knowledge 

of functions An,a for Re(s) G (0,1] implies knowledge of these funetions for s with 
Re(s) > 0. 
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3. Generating functions and integral representation Riemann zeta 

FUNCTIONS AT INTEGER VALUES 

OO 


{ / \ N OO 

> i.e. fn{x) 
^ J j=0 

= We have the following simple observation: 

Proposition 2. Vz £ (—1,1) : f-i(x) = 1, fo(x) = 1/(1 — x) : 

fn(x) = — -^ / fn-l{y)dy, 

Jo 


(3.1) 


c(l — x) 


n > 1. 


ii) Let us denote Q{x,y) the generating funetion of function series {/n}„>o *-e- 
Q{x,y) = XXoV^fjix), fory G (-1,1). l^e have 


(3.2) 


Q{x,y) = 


B{x,l-y,l + y) 


x'^-y{i-xy+y ’ 

where B{x,a,b) denotes incomplete beta function. 

Proof i) We have fn{x) = Xk=i = 

Xk=i XT=k E.=1 x^-^Ht~"^lk = 


■‘3=k 1-a: 

An-l) rs u_i 




dy 


1 V^°° TT(n-l) fS _ 1 rx^oo 

x(l-x) Z^k=l ^k Jo J “y x(l-x) Jo Z^k=l y 

= 7{^fo fn-i(.y)dy. 

ii) We have: {l-x)xQ{x,y) = XT=oy^i^-^)^f 3 i^) = y^ + XT=iy^ Jo 

x+fj X^i y^ fj-i{^)dz) = x+y Jg Q(z, y))dz. Differentiating with respect to x we 
get: (1 — 2x)Q{x, y) + x{l — x)Q'{x, y) = 1 + yQ{x, y). Now solving this differential 
equation we get Q{x,y) = ^eta(xT-yT+yP-C(y) ^ Q(fJ,y) = 1/(1 - y) 

we see C{y) =0. □ 

Let us denote for simplicity C,{s) = X)^i(~l)'’7/* Re(s) > 0. Notice that 
following (1^ we have 


(3.3) 

for fc = 1, 2,... . 
We also have: 


C(fc) = 


r.1/2 


fk-i{x)dx = -/fc(l/2), 


X! = ^{1/2, l-y,l+y), 

3=0 


for y £ (—1,1) following (13.21) . 

Recall that XXiCi‘2j)t^^ = 1 — 7rtcot(7rt) hence XXiCi‘2j)t^^ = ^ 


after some algebra. Hence 

OO 

^ y^^+Xi2j + 1 ) = 5 ( 1 / 2 , l-y,l + y) + l- 

^ sin( 7 rj/) 


sin(7rt) 


- 1 


3=0 
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since C(O) = 1/2. Let us remark that there exist some expansions of incomplete 
beta function. Applying one of them one can we have for example: 


OO 


3=0 


2.-1^ 

3=0 


{-y)j 

j'-U +1 - ’ 


for y £ (0,1). 


4. Remarks on particular values 


In [B] the sums of the form 5(71, fc,/) = E“-oo (jkii)^ > S{n,k,l) = 
were analyzed and some of them were calculated. From the results of this paper it 
follows that the following sums: 

have values of the form tt^ times some know, analytic number. Notice that this 
statement is trivial for k odd, m even and i = 777 / 2 . 

In particular we get for k = 21 we have l/{2m)) — 

({21, {m + i)/{2m))-({2l, {m - i) / {2m)) + ((21, {2m-i)/{2m)) = :^{({2l,l/{2m)) 
+ ({21,1 — 7 /( 2777 )) — ({21, (777 + i)/{2m) — ({21, {m — i)/{2m)). 

Following [B] we also have for fc > 1: 

1 (' 92 fc _ 1 

S{2k, 4,1) = -^{({2k, 1/4) + ({2k, 3/4)) = 2(2fc)! 

where B 2 k denotes 2k — th Bernoulli number. In particular we have 


16K =(C(2,1/4) - C(2, 3/4)); (C(2,1/4) + C(2, 3/4)) = 277^. 

Finally let us recall that ({21, 1) = B 2 i^§^jr■ Using formula (12.51) we get: 


{-lY+^B2i 


(277)2' 

2(20! 


22/-1 
22'-i - 1 


E 


^(2/ —1) 
n2^ 


and consequently we obtain the following expansions of even powers of tt : 


TT 


21 


(-i)'+i 


(20! 

(22'-i - 1)52, 


E 


rr(2^ —1) 

-tin 

n2^ 
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